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O ■ Abstract 
(N 

Neutrino bremsstrahlung of flavor-degenerate pairs in the field of a nucleus 
■ is of potential importance for neutrino astrophysics and is representative of 

a class of processes connecting leptonic electroweak sectors to real or virtual 
photons. We focus on first generation flavor production by both electron and 
\ muon neutrinos and present Standard Model cross sections and distributions 

for lead and iron nuclei. The results (of order 10~^^ cm^ for Vf. - lead collisions 
\Q • at 100 MeV) have been fitted to empirical formulae that can be used to esti- 

. mate backgrounds to neutrino detection experiments and flux normalizations. 

A compact form of the matrix element obtained by analytic reduction is used 
\ to explain the distributions. The V-A limits of the cross sections are shown to 

Qh' agree with published work from the pre-neutral current era. Event signatures 

Q ■ and the possible roles of these processes in stellar and laboratory neutrino 

physics are discussed. Cross sections are compared with those for neutrino- 



> 



C§ . electron scattering for neutrino spectra corresponding to typical supernovae 

temperatures. 

c3 ■ I. INTRODUCTION 

The importance of neutrinos in our fundamental understanding of the physical universe 
is well recognized. Neutrinos, as neutral leptons, interact with matter only via the W and Z 
members of the electroweak gauge boson family. The intrinsic properties of the neutrino play 
a critical role in the evolution and dynamics of the universe and neutrino reactions influence 
the life cycles of stars from infancy to possibly explosive death. Through their dynamic role 
in stellar astrophysics and cosmology, neutrinos also affect the production and distribution 
of elements and serve as valuable informants on the cosmos and its stellar inhabitants. As 
elementary particles and laboratory probes, neutrinos can unravel nature and matter at 
various levels of structure from atoms and nuclei to valence and sea quarks. 

The established importance of neutrino-nucleus and neutrino-lepton electroweak pro- 
cesses in stellar and laboratory physics has motivated increased interest in achieving higher 
accuracy in the theoretical and experimental determination of these rates. This interest 
has stimulated exploration of the possible effects of related neutrino induced electroweak 
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channels that may be competitive to the leading channels or affect their observation and 
detection 

Bremsstrahlung of charged lepton-antilepton pairs by neutrinos in the electromagnetic 
field of nuclei is representative of a class of processes, which produce lepton-antilepton pairs 
in the electroweak sector and couple these electromagnetically to nuclei. As the nucleus 
participates via virtual photon exchange, the resultant cross sections are of order (GZa)^, 
where G is the Fermi constant and (Za) represents the electromagnetic coupling to the 
nucleus. These may be compared to processes of order (G^a), including neutrino-photon 
collisions [^-|^ , and alpha order corrections to the lowest order electroweak channels [0[7|,^ . 
The contributions from pair production could therefore surpass the (G^a) ones for values of 
Z exceeding 12, depending on the other effects such as phase space. They are also compa- 



rable with neutrino-electron scattering channels |]T0| at sufficiently high Z values and could 
be manifested as backgrounds to neutrino cross section and flux measurements as well as 
neutrino detection experiments. The leading order neutrino-nucleus cross sections, behaving 



as (G^A) |TT|-[14| yield larger cross sections than those from pair-production channels. 



The neutrino bremsstrahlung in the nuclear field can produce flavor-asymmetric or sym- 
metric pairs of all flavors, depending on the incident neutrino energy and flavor. Restricting 
exit channels to flavor degenerate charged pairs and including all flavor varieties for the 
incident neutrinos, the allowed channels are: 

+ u'^ + A + 1+ + 
i/^ + A^ z/^ + A + /+ + /-, 

+ ul + A + l^ + 1-, (1) 

where and /~ refer to the produced charged lepton-antilepton pair, and A to the nucleus. 

In this paper, we focus on the production of first-generation flavor pairs only and consider 
incident neutrinos of electron and muon flavor. The tau neutrino rates would be the same as 
those for the muon neutrino. The formalism is in the framework of the Standard Model and 
valid for all combinations of flavor-degenerate and non-degenerate pairs with appropriate 
choice of parameters. Due to the lower threshold, flavor-degenerate first-generation pairs 
are those expected to be of importance for stellar astrophysics or low-energy laboratory 
neutrino physics. The electron neutrinos produce electron-positron pairs by both charged 
and neutral current channels, while the muon neutrino produces these pairs only via neutral 
current interactions. We restrict our study to spin-zero nuclei that respond coherently and 
elastically. 



Processes of this type have been calculated in the pre-neutral current era p!5|-p0|, but 
the literature displays small disagreements between the various cross sections. In light of 
the topical importance of stellar and laboratory neutrino physics, we were motivated to 
provide a correct description of these processes, including the neutral current channels, 
and assess their astrophysical and laboratory importance. We have computed the cross 
sections for these processes in the Standard Model as a function of incident neutrino energy 



and compare the V-A limit with earlier work [^], as the earlier work is restricted to this 
approximation. We have fitted the numerical Standard Model cross sections for production 
of electron-positron pairs by electron and muon neutrinos to determine empirical formulae 
that can be used for various applications. Angular distributions for all final state leptons 
and energy distributions for the produced pair are also reported. The coupled electroweak 
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and photon vertices involved in the processes under investigation are also of interest from 
the point of view of fundamental physics and the parameters of the electroweak theory. 
We use analytic approximations to the cross sections in order to understand these physics 
issues and explain the energy and angular distributions. Possible roles of the reaction in 
stellar and laboratory neutrino physics are discussed. The results can be used to estimate 
effects of the process on neutrino science. The remainder of this paper is divided into the 
following sections: Bremsstrahlung of electron flavor pairs in neutrino-nucleus collisions; the 
formalism; Standard Model cross sections; analytic reduction; distributions; experimental 
signatures and possible roles in neutrino science. 



II. BREMSSTRAHLUNG OF ELECTRON FLAVOR PAIRS IN 
NEUTRINO-NUCLEUS COLLISIONS 



Electron-flavor pair production in neutrino-nucleus collisions can proceed by both 
charged and neutral electroweak currents. The incident neutrino flavors allowed by the 
Standard Model are necessarily different for the two types of currents, and the Feynman 
diagrams for such processes include a purely leptonic electroweak sector coupled to the nu- 
cleus via a virtual photon. The channels can be represented by the general equation Eq. (1). 
The Feynman diagram for charged current mediated electron-positron pair bremsstrahlung 
is shown in Fig. 1 along with the cross term for electron neutrino-initiated channels. The 
neutral current processes allow flavor freedom to the incident neutrino and can be described 
by the Feynman diagram of Fig. 2. It may be noted that the electroweak sectors in Figs. 1 
and 2 represent purely leptonic currents interacting via the electroweak gauge bosons. In 
the absence of the nuclear legs and their intermediary virtual photon, the electroweak sector 
would describe the process, i/ — z/ + e~ -l- e"*", analogous to muon decay. Such a channel 
would satisfy Standard Model charge and lepton-number conservation by a suitable choice 
of the participating lepton flavors and charges and could be compared with other purely 
leptonic four-fermion electroweak channels. However, it cannot be realized physically due 
to kinematic and momentum constraints. 

The presence of the virtual photon and nuclear sector make the pair production process 
on-shell, and its interference with the electroweak sector determines the correlations and 
physics of these processes. We demonstrate the presence and survival of the individual 
sectors via an analytic reduction of the matrix element. This also permits analysis of the 
distributions in the light of helicity constraints similar to those appearing in purely leptonic 
electroweak modes. The cutting of the nuclear legs in Figs. 1 and 2 transforms the virtual 
photon into a real one, and the diagrams then correspond to similar pair production channels 
in neutrino-photon collisions, which can be described by an analogous formalism. We discuss 
these issues further in a later section. In contrast to the processes described by Figs. 1 and 
2 and Eq. (1), the neutrino-nucleus reactions believed to dominate explosive star death are 
the lowest order electroweak processes corresponding to direct exchange of W or Z gauge 
bosons between the neutrino and the nucleus and contributing to order (G^A) 



III. THE FORMALISM 
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A. The Electroweak Sector and the Effective Lagrangian 



Neutrinos produced in supernova environments, solar interiors, and those involved in 
laboratory experiments have relatively low energies, and this restricts the four-momentum 
transferred at the electroweak sector to values well below the masses of the electroweak 
gauge bosons. Under these circumstances, the gauge bosons may be collapsed to point 
coupling. As we are primarily interested in such neutrinos, we will carry out the numerical 
calculations in the effective Lagrangian approximation of the Standard Model using the 
reduction of the gauge boson propagators to local couplings. Since the electroweak vertices 
in Figs. 1 and 2 are purely leptonic, the reduction to the low-energy, effective Lagrangian 
coupling is straightforward. 

The heavy photon-like W-boson propagator of mass Mw can be approximated by []TU 
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W 



when the four-momentum transfer qy/ is such that >>| |. The dimensionless 

electroweak coupling constant g of the Weinberg theory can be related to the universal 
Fermi constant G through 

G/v^ = W(8M2.), (3) 

and the effective Lagrangian for neutrino-electron coupling via charged electroweak currents 
can be expressed as 

Llff = ^m^m^ - 75)^(^)][V5(e)7,(l - 75)^(e)]. (4) 

The corresponding effective local second-order S-matrix for neutral current-induced elec- 
troweak processes has the form: 

L:^f = ^mumi - ^,)^{umeh,{a - h^me)] (5) 

where 

a = (—1/2) -|- 2sin'^9w, b = (—1/2), and 9w is the weak mixing angle [^,0. Channels 



proceeding via both charged and neutral current corresponding to the z/g-induced {e~^e~) 
pairs, can be described by a combined effective Lagrangian of the type Eq.(5) given above, 
with a ^ a + 1 and b ^ b + 1. -induced (e"'"e~) pair production is described by Eq. (5). 
The V-A limit can be reached for a ^ 1 and 6—^1. These limits are given in Table |. 

Using collapsed propagators, the electroweak sector of the diagrams in Figs. 1 and 2 
reduce to point coupling, and the charged and neutral current channels acquire apparent 
degeneracy as shown in Fig. 3. Figure 3 also describes the (V-A) limit. We note that the 
Fierz rotation has been applied to the charged current channel as customary. 
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B. The Nuclear Sector and Its Coupling 



The nucleus participates via virtual photon exchange with one of the members of the 
produced flavor-degenerate pair. The virtual photon transfers four-momentum q between 
the lepton propagator and the nucleus, where q is defined through q = p — p', and p and p' 
are the initial and final four-momenta of the nucleus. Since the nucleus is connected to the 
electroweak vertex through the virtual photon, the nuclear interaction is purely electromag- 
netic and only the electromagnetic form factor of the nucleus is operative. As mentioned 
earlier, we only consider spinless nuclei and collisions that leave the nucleus in it's ground 
state. The nuclear matrix element of the electromagnetic current for a spinless nucleus can 
be expressed as 

<p'\jLiO)\P>=iP + P'f ZeF{q'). (6) 

The form factor limits the four-momentum transferred to the nucleus to be less than the 
approximate reciprocal size of the nucleus, q < 20 MeV/c. It has been shown by that 
allowing q its full range of allowed kinematical domain yields an erroneous E^/n E asymptotic 
dependence to the cross section instead of the correct E /n E dependence obtained using the 
form factor description. Our studies validate this observation. We use an exponential form 
factor given by 



F{q')=exp[-\q'\/4), (7) 

where, kq = 116.46^1 MeV/c. 

While we have confined this study to cases where the nuclear response is coherent and 
elastic, the nucleus in all of these processes could respond inelastically. As energies in- 
crease, spontaneous release of nucleons or nuclear excitations accompanying the neutrino 
bremsstrahlung by pair radiation may contribute to the total cross sections. We do not con- 
sider incoherent contributions from individual nucleons in the present work. At the energies 
we are interested in (1 to 1000 MeV), these are expected to be smaller than the elastic cross 



sections. Contributions from incoherent processes have been studied in earlier work |15,20 



C. The Matrix Elements and Cross Sections 

The second-order S-matrix elements corresponding to electron-positron pair production 
by neutrinos can be written in a generalized form in terms of the appropriate effective 
Lagrangian as 

M = -^iZe')Fiq')^[Uik'mi - l.Mk)] 

where 

= — q, = g — r+ and Q refers to 7aQ" in general. Values of a and b have 
been defined earlier and differ for the separate cases of incident z/^ and z/g and for the V-A 
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case. Equation (8) includes the direct and cross terms. U and V refer to the spinors for the 
respective particles or antiparticles, and and (5+ are the four-momenta of the charged 
fermion propagators. The four-momentum transfer to the nucleus is g, and e is the usual 
electromagnetic coupling constant. P and q are connected to the four-momenta of the initial 
and final nuclei through q = p — p' and P = p + p'. F{q^) represents the nuclear form factor 
as discussed earlier and controls the momentum transferred to the nucleus at high energies. 
It is unity at low energies. The above labeling of the momenta is shown in Figs. 2 and 3. 
Introducing the phase space factors, the cross section can be written in terms of the matrix 
element as 



a 



+ ^ \M\^ j^^5\k' + r+ + r_-k-q) (9) 



AE,E,J 2E'2E+2E.2E,,^f-^' ' (27r) 



Conventions for phase space factors and gamma matrices are those of [^. The square 
of the matrix element is summed over final spins. There is no averaging over initial spins, 
as the neutrino is helical. The delta function ensures four-momentum conservation for the 
process. 

Including flavor variety, the formalism can be readily extended to degenerate flavor pro- 
duction beyond the first generation and to nondegenerate flavor production using appropriate 
choices for the constants a and b. 



IV. STANDARD MODEL CROSS SECTIONS 

The three final leptons and the recoiling nucleus span a twelve dimensional final state 
phase space. This twelve-dimensional phase space is reduced to eight- dimensions by the 



delta function constraints. The phase space reduction is carried out in the usual way 
considering the charged pair in their center of mass frame first and then incorporating 
the final neutrino. One azimuthal angular integration is removed by symmetry arguments 
and a suitable choice of axis. The remaining seven-dimensional phase space is evaluated 
numerically using Monte Carlo techniques. 

We display in Fig. 4 a comparison between the earlier results for this cross section 
obtained by Cyz et. al. |15| and the charged current results for i/g-induced events in nickel. 



The latter corresponds to the V-A approximation used by [|I5| and the nuclear form factor 
used therein to enable the comparison. For our Standard Model results presented in all 
other figures, we use the form factor discussed in the previous section. 

In reality, the correct cross sections should include neutral current contributions. Figure 5 
displays the cross sections (including neutral currents) for and induced channels for ^^Fe 
and ^°^Pb. The points represent the values of the numerical computations corresponding to 
up to 10^ Monte Carlo points and a statistical error better than five percent. The induced 
cross sections are lower than the z/g initiated ones as expected because the former channels 
are mediated only by neutral currents, while the latter include both charged and neutral 
currents. Clearly, the contributions of both processes increase with Z and are expected to 
be most important for high-Z nuclei. Iron and lead have been selected as representative 
nuclei for medium and high Z nuclei that are important in supernova neutrino physics and 
laboratory experiments, including the terrestrial detection of cosmic neutrinos. A naive 
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scaling of the results for iron could be obtained to obtain cross sections on lighter nuclei like 
oxygen or carbon. 

We have also fitted the energy and Z'^ dependence of the cross sections so as to represent 
them by empirical expressions that can be used to facilitate applications and estimate the 
effects of these channels in different environments. The numerical expressions for Vf. and z/^ 
induced channels are given below and were obtained as an empirical fit to the numerically 
computed cross sections in iron. 

The cross sections, in units of cm^, for first generation pair production by muon- and 
electron-neutrino bremsstrahlung in iron can be represented as, 

ln((T) = C^x ln(x) + C2X + C3 (10) 

where x = \nE^ and E,^ is the incident neutrino energy, in MeV. 

The coefficients Ci, C2, C3 are given in Table ||, and comparisons with numerical calcu- 
lations are shown in Fig. 5. The solid lines in Figure 5a correspond to values given by the 
fitted expressions in Eq.(lO) with values from Table |I[ A scaling of the above expression 
gives excellent fits for both carbon and lead, demonstrating their usefulness for various ap- 
plications. The solid curves in Fig. 5b display values obtained using this expression, scaled 
from iron to lead. 



V. ANALYTIC REDUCTION 

The complexity of the expression for the spin summed square of the matrix element, 
with its many terms does not lend itself to analysis of the energy and angular distributions 
or to the interpretation of the interesting physics of such combined electroweak and photon 
exchange diagrams. These features however acquire clarity when examined via a reduced 
form of the matrix element. 

The matrix element of Eq. (9) can be reduced to a compact form by using the Dirac equa- 
tion. Rationalizing the fermion propagators and replacing and (5+ by their expressions 
in terms of r_ and r_|_, equation (8) takes the form 

M = ^{Ze')F{q')^[U{k')r{l - l,)U{k)] 

- {r4-4 + m) u \ ^ ( h 11./ X /I IN 

^ (^~)[^" {^Ql-rn?) ~ ^^5) + 7m(« - ^75) (g2__^2) 7jV^(r+) (11) 



Rearranging, using usual spinor algebra , operating on the electron and positron four 
momentum vectors from the left and right by their respective spinors, and using the Dirac 
equation, the matrix element can be expressed as the sum of two terms M = Mq + M', 
where 

Mo = ^^[t/(fc')7^(l - 75)f/(fc)] 



X 2 



V2 q' 

(Pr_) (Pr^ 



{Ql - m2) {Ql - m 



2} 



[f/(r_)7^(a-&75)V^(r+)] 

(12) 
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and 



M' 



X f7(r_)[- 
X 1/(r+). 



-f/(fc')[7''(l-75)f/(fc)] 
^7/^(0 - ^75) 7m (a - ^75 



(13) 



This decomposition allows us to isolate the part that coincides with the standard electroweak 
leptonic processes like v — e scattering and muon decay, for which the spin summations and 
angular correlations are well known. The nuclear sector with its virtual photon connector is 
also identifiable. Mq allows us to attempt an analytic understanding of the physics of such 
interfering vertices, which would not otherwise be realizable due to the underlying numerical 
complexity of the terms. The remaining term M' is not amenable to such a compact form. 
Before assessing the contribution of M', we study the behavior of Mq. Neglecting M', the 
cross section takes the form: 



S'k' (fr^ d?r^ S'p' 



AE„E^ J 2E' 2E^ 2E^ 



2Er,l 

P spins 



E 1^0 



{2n) 



The spin summed value of | Mq p has the form 



E 1^0 



spins 



2\2 



{Pr, 



(14) 



(Q 



x64[(a + hf{kr+){k'r^) + (a - hf{kr^){k'r+) 

W(a2 -6')(A;'A;)] 



(15) 



The result of the trace part of this is readily seen to be identical to pure leptonic four- Fermi 
processes and can be compared by inspection. (The sign of the mass term differs from that 
in neutrino-electron scattering due to the occurrence of a final state positron instead of an 
incident electron.) 

We present in Fig. 6, a comparison of the exact numerical calculations with the results 
obtained using the Mq term from the analytic reduction. The agreement indicates the 
neglected terms do not contribute appreciably, and thus the simplified expression of Eq. 
(15) may be used for various applications. Expression (15) accounts for about 80 percent 
of the total cross section. The connection between the processes studied in this paper 
and photo-neutrino pair production can be easily seen by cutting the nuclear legs in the 
Feynman diagrams of Fig.s 1-3 |25]. In this case, there is no equivalent of the term M', so 
that one could conclude that in the limit — 0, the term M' should vanish. The fact that 
usually virtual photon processes can be well approximated by photo-nuclear cross sections 
or the equivalent photon method p2| , p!6| , p!9| vindicates our numerical demonstration of the 
dominance of Mq in M. While all cross sections reported here contain the complete matrix 
element, we have used the reduced ones to compute and explain the distributions reported 
in the next section. 

The survival of the leptonic electroweak sector in factorable form in Mq allows us to 
compare the pair bremsstrahlung with the expressions for neutrino-electron scattering and 
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other lowest-order electroweak channels. The angular constraints and distributions typical 
of weak processes can also be recognized through this sector. The factorable nuclear sec- 
tor imposes the physics constraints typical of virtual photon-exchanged processes, and the 
ultimate distributions and energy behavior are determined by the interference of the two 
sectors. 



VI. THE DISTRIBUTIONS 

Angular distributions for all the emitted leptons are shown in Fig. 7 for electron neutrino- 
lead collisions at 10, 100 and 1000 MeV incident energies. The angular distributions are seen 
to peak in the forward direction, the outgoing neutrino having the sharpest rise. All the final 
leptons display a sharpening of the spectra, emphasizing their forward bias with increasing 
energy and the angular distribution of the positron is observed to be broader and less forward 
peaked than that of the electron. The characteristics are discussed below in the light of the 
conflicting dictates of the electromagnetic and weak vertices. The angular distributions for 
this process are reported here for the first time and can be used, along with the energy 
distributions and cross sections to estimate experimental signals. 

Energy distributions of the charged leptons are displayed in Figs. 8 and 9 for lead nuclei. 
The general nature of these distributions can be explained in terms of the high-energy limit of 
the analytically reduced matrix elements and can be shown to exhibit the helicity constraints 
of the electroweak vertex. 

Figure 8 provides a comparison of electron and positron energy distributions for incident 
neutrinos of 10, 100 and 1000 MeV energies. The horizontal axes are scaled by the incident 
energy. In general, the energy distributions display a differentiation in the spectra of the 
two members of the produced charged lepton pair, the positron spectra peaking at lower 
energies and more sharply than that of the electron. The effect is increasingly pronounced 
with increase of incident and available energy, the positron peak sharpening and shifting to 
lower energies. 

Figure 9 compares the electron and positron distributions for 100 MeV incident neutrinos 
of muon and electron type respectively with those for the V-A case. It may be noted that the 
differentiation in the lepton-antilepton spectra displayed in Fig. 8 for the electron neutrino 
events is retained in the data of Fig. 9, although the effect is more pronounced for the 
V-A case and least for the muon neutrino channels. The asymmetry between the electron 
and positron spectra obtained by us is consistent with a similar asymmetry reported for 



the V-A case by |^0[ for second generation pairs produced at higher incident energies. The 
contributions of the neutral current sectors soften the differentiation of the spectra of the two 
charged leptons for the Standard Model distributions in electron neutrino induced events 
in contrast to the V-A case. This is why the asymmetry is most pronounced in the V- 
A limit, corresponding to charged current events only. The muon neutrino produces first 
generation pairs by neutral current only and the differentiation in the energy distributions 
of the outgoing pair is barely perceptible. 

To understand the distributions let us first look at the high energy limit of the spin 
summed value of | Mq p, in particular the various four vector products. We can write 
I Mq P from section V, equation (15/14) as 
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where 



and 



I Mo |2= 



spins 



2 (g^)^ 



(16) 



(17) 



D2 = 64[Ai{kr+){k'r_) + A2{kr^){k'r+) + A^{k'l 



(18) 



with Ai = {a + 6)2, = (a - 6)^ and A3 = (a^ - 1)^)711? 

The constants a and b have the values stated in earher sections for specific neutrino 
flavors, and in Table |. 

We examine the first term, containing four vector products (A;r+) and {k'rJ) multiplied by 
the coefficient Ai. In the high-energy limit, for relativistic pairs, the first four- vector product 
(fcr+) of this term approaches (fcr+) — >• £'^£'+[1 — cos(^+) -|- Oiw?)] where 6*+ represents the 
angle of the produced positron with respect to the forward direction or the direction of the 
incident neutrino. This term goes to zero as cos (6'+) goes to unity, and therefore suppresses 
relativistic positrons in the forward direction. The product refiects the helicity transfer at 
the weak vertex and originates in the charged current component of the combined process. 
The behaviour is reminiscent of similar helicity constraints that infiuence neutrino- and 
antineutrino- electron scattering at high energies and distributions in muon and pion decay 

El- 

The suppression of relativistic forward positrons, dictated by the weak vertex is opposed 
by the virtual photon sector that links the weak vertex to the recoiling nucleus. The presence 
of in the denominator of Di yields the largest probability for those events with the 
smallest values of q^. This tends to focus all final leptons into a narrow forward cone at high 
energies, since energy-momentum conservation takes the form: kK.k' + r^ + r_^ for small 
q^. This is common to many other processes at high energies, including neutrino-electron 
scattering and pair production [^. As q^ occurs as a square in the denominator it is 
powerful in maintaining the clustering of all final leptons into a narrow forward cone, with 
an opening angle of the order of {m/EyY. The analytic expression for the spin-summed 
square of the matrix element allows us to understand the competing effects that infiuence 
the final distribution. The confiict between the forward cone dictates of usual high-energy 
exit channel behavior and the strong suppression of forward positrons by the weak vertex 
and its helicity constraints can be best compromised by suppressing high energy positrons 
and enhancing positron population of lower energy states. This behaviour is similar to 
the suppression of pion decay by the electron mode despite its favoring by phase space as 
compared to the dominant muon mode. 

The above arguments are supported by the angular and energy and distributions in 
Figs. 7,8, and 9. All final leptons dominantly go forward in accordance with traditional pair 
production observations, and this trend is enhanced at higher energies. On the other hand, 
to satisfy the helicity requirements of the charged current component, the positron energy 
distribution peaks at lower energies and is strongly suppressed at high energies. This is the 
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direct effect of the term kr^ which goes to zero for relativistic forward positrons. As incident 
energies increase, the charged leptons would normally tend to become increasingly relativistic 
and the differentiation in their spectra due to suppression of relativistic forward positrons 
is accentuated further as seen in Fig. 7. This behavior is also consistent with helicity 
arguments similar to those that explain the asymmetry in neutrino-antineutrino scattering 
at high energy, correlations in muon decay and other weak interaction processes. The second 
terms do not display this requirement and their interference softens the asymmetry in the 
distributions as compared to the V-A case. 

Turning to the electron spectra, the second four vector product {k'r_), involving the 
electron and the final neutrino, reflects the helicity requirements on these leptons and a 
similar suppression on their collinearity. However, this is different from the previous case, 
where the incident neutrino defines the forward direction. For the electron-final neutrino 
correlation term, both electron and final neutrino have freedom of direction. This allows 
a spread to the correlation between them. However, this too suppresses forward electrons, 
collinear with the outgoing strongly forward neutrino, and hence peaks the electrons at 
intermediate energies. The forward suppression is milder for the electrons as it occurs 
relative to the final neutrino. This differs from the case of the positron, which is suppressed 
with respect to colinearity with the incident neutrino which in fact defines the forward 
direction. This can explain the sharper and lower energy peak for the positron as compared 
to the electron and the stronger forward production of electrons. 

We now investigate the importance of the first term that favors lower positron energies 
compared to the others in D2 and through it on the expression for a by quantitative compar- 
isons of the coefficients. Table |T| displays the numerical values of the coefficients Ai and A2 
in the expression for D2 for the V-A, u^. and z/^ cases. (The interference term, characterized 
by the coefficient is mass dependent and does not contribute much at high energies). 

It is apparent from Table |T|, that Ai dominates for Ue and pure charged current (or 
V-A) events, the dominance being more for the V-A case when the second coefficient is 
zero. This explains the strong differentiation in the charged lepton spectra for these classes 
of events. The presence of the second term, arising from the interference from the neutral 
current contributions, softens the z/g spectra slightly, compared to the V-A case. On the 
other hand, for i/^ events, both coefficients are of similar magnitude, giving rise to the much 
softer differentiation in the charged lepton spectra for these events. 

The energy and angular distributions reported here have been calculated in the labo- 
ratory frame using the reduced matrix elements to enable comparison with corresponding 
analytic expressions. A check with the complete matrix elements, indicated that the agree- 
ment between the reduced and complete matrix elements is even better for the distributions 
than for the cross sections. Further, the distributions are of use and interest mainly to dis- 
play the underlying physics and to be indicative of experimental signals. This contrasts our 
interest in the cross-sections, which are useful to obtain numerical estimates of the influence 
of these processes on neutrino physics in general. We reiterate that cross section results in 
Figs. 4-6 include the full matrix element as do the empirical fits that represent their energy 
dependence. 

VII. EXPERIMENTAL SIGNATURES AND POSSIBLE ROLES 
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IN NEUTRINO PHYSICS 



From the laboratory perspective, the event signature would be the produced pair, both 
members being emitted predominantly in the forward direction. The opening angle of the 
outgoing positron is larger than that of the electron as discussed earlier. Further, the 
positron energies peak at the lower energy end of the spectrum while electron energies have 
a broader distribution of energies. Events of this type could be studied by detecting both 
components of the produced pair in coincidence. Such pairs should be distinct from those 
arising from de-excitation of nuclei, as the energy available to the latter pair is fixed by the 
de-excitation energy of the nucleus and the correlation signature of such a de-excitation pair. 
In contrast, the energy distributions of the neutrino- initiated direct pairs display continuum 
characteristics, and the angular distributions respect helicity constraints of the electroweak 
vertex and the physics of the photon vertex. 

If the primary experimental interest is not to detect the bremsstrahlung process itself, but 
rather its influence on measurements of the lowest order neutrino-nucleus events or neutrino 
flux determinations, this can be estimated using the theoretical cross sections presented in 
this paper and the empirical formula obtained for them. A major focus of current neutrino 
science is the detection of neutrinos of extra-terrestrial origin, - those from the sun or from 
supernovae. The primary event signatures for these are from charged current neutrino 
interactions occurring on the detector target nuclei. Detection processes that measure the 
outgoing electron in neutrino induced charged current events could pick up pseudo-events 
generated by electrons from bremsstrahlung events, particularly in high Z targets like lead 
or iodine. As experimental constraints tighten and increased accuracy in event rates is 
attempted, it becomes important to estimate the effects of backgrounds such as the pair 
bremsstrahlung investigated in this paper. As mentioned earlier, the pair bremsstrahlung 
cross sections of order {GZaY may become comparable to order alpha corrections to the 
lowest order electroweak channels and G'^a order processes for values of Z exceeding 12. 

The process could also be a potential source of background to neutrino-electron scattering 
experiments and experiments that use u — e scattering to establish flux measurements and 
normalizations. This background will be increasingly important for higher Z values of the 
nuclei of the material hosting the participating electrons. Even for light nuclei, the process 
could be a source of background for experiments that seek fiux-accuracy to a few percent, 
or better. 

We present pair bremsstrahlung cross sections folded with neutrino energy spectra for 
neutrinos emanating from thermal sources with typical supernovae temperatures in Table 
IV| . We compare our results for iron and lead with those for neutrino-electron scattering. We 
have used for the neutrino distributions, the normalized Fermi-Dirac form given in |T2| for 



appropriate supernova temperatures. As expected, the contribution of pair bremsstrahlung 
increases as a relative percentage of the neutrino-electron scattering for increasing Z and 
energy values. 

The results in Table |TV| include thermal sources at temperatures higher than those ex- 
pected to generate the bulk of electron neutrinos, which decouple at temperatures lower 
than the higher flavor neutrinos. This inclusion allows for their use in applications involv- 
ing supernova thermometry that explore oscillations from tau and muon to electron flavors 



12,11 . Muon or tau neutrinos that have oscillated to electron flavor would interact in a 
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terrestrial detector according to the electron neutrino values at the appropriately higher 
temperatures. 

A comment about anti-neutrino events is in order. The diagrams of figures 1-3 can be 
easily modified to describe anti-neutrino induced events. As a particle-antiparticle pair is 
produced, crossing symmetry dictates identical cross sections for incident neutrinos or anti 
neutrinos. This has been discussed in additional detail in [jl5| . However, the distributions in 
figures 7-9 would display an interchange between the charged leptons due to the well-known 
fiipping of the charged current part (appearing through the first term in D2) as neutrino- 
electron scattering changes from incident neutrinos to antineutrinos. This does not alter 
the final cross sections and the results presented here for neutrinos would apply equally to 
antineutrinos. 

Apart from the infiuence of pair bremsstrahlung on the detection of cosmic neutrinos and 
hence on our understanding of the astrophysical events they portray, the question of possible 
roles in the astrophysical dynamics themselves also arises. What might be the astrophysical 
environments to spawn such neutrino-bremsstrahlung pairs? High Z environments such 
as pre-collapse supernova cores, scenes of the stalled shock, or the neutrino-driven winds 
hosting the r-process and neutrino nucleosynthesis suggest themselves as potential sites. 
The density of the Fermi sea, however, would inhibit low-energy electron production in pre- 
collapse cores. It may be noted that the Coulomb distortion effect would affect low-energy 
electron and positron emission oppositely and could be considered, along with the Fermi sea 
effects, to determine the cross sections in physical supernova environments. 

The question does arise whether the pair bremsstrahlung might be more important for 
anti neutrinos in supernova environments with a high population of neutron-rich nuclei. As 
the lowest order nuclear channels are suppressed due to non- availability of neutron states 
as well as the negative effects of the Coulomb barrier |TT[|, other channels such as pair 



bremsstrahlung may provide alternate means of slowing anti-neutrinos. 
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TABLES 

TABLE I. Given in the table are the dimensionless parameters, a,b, of the effective Lagrangian 
for different flavors of the incident neutrino. 



Neutrino Flavor 


a 


b 




2sm^{Qw) - \ 


1 

2 




2sin2(eiy) + i 


+ ^ 


V-A 


+1 


+ 1 



TABLE IL Values of the coefficients Ci, C2, C3 used in Eq.(lO) for representing the total pair 
Bremsstrahlung cross sections for both electron and muon incident neutrinos on iron. 



Neutrino Flavor 


Ci 


C2 


C3 




0.34032 


5.15593 


109.5118 




0.37696 


4.64135 


110.4906 


TABLE III. Given in the table are the dimensionless parameters, A\, and A2 of (Eq.s (16) - 
(18))for different flavors of the incident neutrino. 


Neutrino Flavor 


Ax 






2.13 


0.21 




0.29 


0.21 


V-A 


4 






TABLE IV. Neutrino induced pair cross sections off of iron and lead for normalized 
Fermi-Dirac spectra of temperature T are compared to those for (fg, e~) elastic scattering. 



T (McV) 


(Te X 10 ^^cm^ 


(7'pair(i^e)/cre 


<7pair(i^b)/o-e 


4.0 


0.113 


0.007 


0.103 


6.0 


0.170 


0.018 


0.180 


8.0 


0.323 


0.026 


0.258 


10.0 


0.328 


0.034 


0.335 


12.0 


0.342 


0.041 


0.412 


15.0 


0.425 


0.053 


0.525 
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FIG. 4. Comparison of neutrino induced pair cross sections, as a function of laboratory neutrio 
energy, off of nickel using the V-A approximation with the earlier results of Ref. 15. 
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FIG. 5. a) Electron-neutrino (solid circles, and solid line), and muon-neutrino (solid squares 
and dashed line), initiated Standard Model pair cross sections on iron evaluated numerically (circles 
and squares) and compared to the empirical cross section formulas, Eq.(lO) (solid and dashed lines), 
b) same as in a) for the nucleus lead, for this case, the solid and dashed lines represent Eq.(lO) 
with an assumed Z"^ scaling. 
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FIG. 6. Electron-neutrino induced pair cross sections off of lead as a function of the labo- 
ratory neutrino energy, computed numerically, and compared those obtained using the analytic 
approximation discussed in the text 



22 



50 
40 -_ 
30 
20 
10 - 
150 



X 
■D 



125 -_ 
100 -I 
75 -I 
50 -I 
25 ^ 


250 ^ 
200 ^ 
150 H 
100 
50 - 




-1 — I — I — r 



T — I — I — r 



+ 



.4 



,o 



I I I I I I I I 



I I I I I I I I I I I I I I 



■— 


-10 MeV 


- o 


-100 MeV 




- 1000 MeV 



P- 



i^rrr i i i i I i i i i 



. Oj 

. / 

*/ 

Jo 



^y-B— ■ 



0.95 



0.96 



0.97 



0.98 



0.99 



1.00 



X 



FIG. 7. Distributions for the electron, positron and the outgoing neutrino, as a function of 
X, the cosine of the outgoing angle with respect to the ingoing neutrino direction, for electron 
neutrinos of incident energies 10, 100, and 1000 MeV on lead. Note that all three distributions are 
sharply peaked in the forward direction, the positron distribution being less so than the others. 
The distributions are normalized so that / dxda/dx = 1. 
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FIG. 8. Energy distributions for electrons and positrons for electron neutrinos incident on lead 
with energies of 10, 100, and 1000 MeV. The distributions are normalized so that / dyda/dy = 1, 
y = E/Ey, where Ey is the energy of the ingoing neutrino. 
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FIG. 9. Energy distributions for electrons and positrons for neutrino incident on lead at an 
energy of 100 MeV for electron neutrinos, muon neutrinos, and for the case (V-A). The distributions 
are normalized so that J dyda/dy = 1, y = E/E^, where Ei, is the energy of the ingoing neutrino. 
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